Crystal-field ground state of the orthorhombic Kondo insulator CeRu 2 Ali 
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We have succeeded in establishing the crystal- field ground state of CeRu2Alio, an orthorhombic 
intermetallic compound recently identified as a Kondo insulator. Using polarization dependent soft 
x-ray absorption spectroscopy at the Ce M^^ edges, together with input from inelastic neutron 
and magnetic susceptibility experiments, we were able to determine unambiguously the orbital 
occupation of the 4/ shell and to explain quantitatively both the measured magnetic moment along 
the easy a axis and the small ordered moment along the c-axis. The results provide not only a 
platform for a realistic modeling of the spin and charge gap of CeRu2Alio, but demonstrate also 
the potential of soft x-ray absorption spectroscopy to obtain information not easily accessible by 
neutron techniques for the study of Kondo insulators in general. 

PACS numbers: 71.27.+a, 75.10.Dg, 75.30.Cr, 78.70.Dm 



CeRu 2 Al 10 is a fairly new synthesized orthorhom- 
bic intermetallic compound 1 ' 2 and has initiated a flurry 
of experimental and theoretical research activities since 
the discovery of its Kondo insulating properties in 
20094 — iiir— The electrical resistivity shows a thermally 
activated behavior at elevated temperatures and has a 
more metallic like form several degrees below To = 27 K 
where a phase transition occurs as revealed by specific 
heat measurements j 3 -^ A large anisotropy is observed in 
the static susceptibility (% >Xc >Xb) and along the easy 
axis the susceptibility is Curie- Weiss like above T .— 
For some time it has been a matter of debate whether the 
phase transition is magnetic in origin^r— i 12 i 13 In fact, in 
the structure of CeRu2Ali (space group Cmcm) each 
Ce atom is situated in a polyhedral cage consisting of 
16 Al and 4 Ru atoms, giving rise to very large Ce-Ce 
distances of more than 5A<2 Thus, it is not obvious to 
explain the relative high ordering temperature using a 
standard model based on RKKY exchange interactions. 

Nevertheless, recent muon spin relaxation and neutron 
diffraction experiments confirmed the magnetic nature 
of the phase transition. The Ce moments order anti- 
ferromagnetically along the c direction with a small or- 
dered moment of [i ord 0.4 ii B Zr&±2. The latest NMR 
data agree with this scenario. 13 Below To inelastic neu- 
tron scattering (INS) experiments on a polycrystalline 
sample have found the existence of a spin gap of several 
meV, characterized by a magnetic mode with flat disper- 
sion at about 8 meVi^ Above To this excitation is drasti- 
cally suppressed and therefore cannot be attributed to a 
crystal electric field (CEF) excitation. More recent INS 
results have found CEF splittings of 30 and 46 meV It 
has been suggested that the mixture of Kondo and CEF 
effects may be responsible for the small size of the or- 
dered moment. A point charge model yields a moment of 



only 0.62 /jb along c for the ground state (GS)jil but the 
combination of large spin gap, small ordered moment and 
high ordering temperature remains quite puzzling. Re- 
cent optical studies have suggested an electronic struc- 
ture with strong anisotropic hybridization between the 
4/ and conduction electrons which is weakest along the 
crystallographic b direction. It has been then speculated 
that a charge density wave can form along b which in 
turn may trigger the magnetic transition ! i 15 

The objective of the present work is to elucidate the 
local electronic structure of the 4/ electrons. The GS 
wave function of the 4/ electrons in this orthorhombic 
compound is expected to be highly anisotropic due to 
the presence of the CEF, with important consequences 
for the magnetic properties and the gap formation as 
pointed out already by several groups in their study 
on other orthorhombic semiconductors like CeNiSn and 
CeRhSb^ 3 --^ 

Efforts to prove these theories experimentally turned 
out to be challenging because broadening of the CEF ex- 
citations due to hybridization and the existence of spin 
gaps prevented the determination of the ground states 
in these compounds by inelastic neutron scattering 22r2L 
Hence there is need for an alternative approach. Our 
method of choice is x-ray absorption spectroscopy (XAS) 
at the Ce M^s edge which has shown to be a comple- 
mentary technique to neutron scattering to determine 
the GS wave function of CEF split tetragonal cerium 
compounds^ - — The sensitivity to the initial state sym- 
metry is achieved via the different absorption for light 
polarized E || c (with c being the long tetragonal axis) and 
E _L c. We now extend this method to an orthorhombic 
compound where the polarization dependence of the ab- 
sorption, the linear dichroism (LD), has to be measured 
for all three directions, i.e. for E \\ a, E \\ b, and E || c in 
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order to determine the initial state symmetry. 

Single crystals of CeRv^Alio were grown by an Al self- 
flux method^ and their quality and orientation were con- 
firmed by Laue x-ray diffraction. All XAS measurements 
were carried out at the Dragon bending magnet beamline 
BL11A1 of the NSRRC in Taiwan. The energy resolu- 
tion at the Ce M 4j5 edge (hv « 870- 910 eV) was 0.4 eV. 
The single crystals were cleaved in situ in an ultra high 
vacuum of ~ 10~ 10 mbar to obtain clean sample surfaces. 
We measured the total electron yield (TEY) and normal- 
ized the signal to the incoming flux Jo as measured on an 
Au-mesh at the entrance of the experimental chamber. 
Two crystals were investigated: one mounted with the c 
and the other one with the a axis parallel to the Poynting 
vector of the incoming light. This way we were able to 
vary the electric field from E \\ a to E \ \ b and from E \\ b to 
E || c, respectively, by rotating the crystals in steps of 90° 
around the Poynting vector. The spectra for E\\b were 
recorded on both samples so that the three polarizations 
could be normalized to each other. The data were repro- 
duced by probing several positions on the samples and 
recleaving the crystals. 

Ionic full multiplet calculations were performed with 
the XTLS 8.3 program^! to simulate the XAS spectra. 
Initially, the experimental isotropic spectra /isotropic = 
lE\\ a + lE\\b + Ie\\c were fitted. The best agreement was 
achieved with a reduction of the atomic Hartree-Fock 
values of about 40% for 4/ — 4/ Coulomb interactions 
and of about 20% for the 3d — 4/ interactions. The re- 
duction factors account for configuration interaction ef- 
fects not included in the Hartree-Fock scheme. Their size 
is in agreement with previous findings on tetragonal Ce 
systems i^2r— We follow the coherent approach for our 
simulations, i.e. the XAS spectra are calculated directly 
from the CEF mixed GS wave functions, the latter being 
fabricated via CEF parameters since for orthorhombic 
symmetry the XAS spectra cannot be calculated as inco- 
herent sums of the pure J z spectra. 

The CEF parameters are defined via the CEF Hamil- 
tonian which arises from the expansion of the Madelung 
potential in spherical harmonics: 

oo k 
k—0 m— — k 



q™(6»,$) = \l 2hTTYl n (6^) are the semi-normalized 

spherical harmonics. The expectation value (r fc ) of the 
radial part of the wave function cannot be calculated ana- 
lytically and is therefore included in the CEF parameters 
AV? — A™(r k ), which are determined experimentally. 

There are nine independent CEF parameters, A®, A|, 
Al, Al, Al, A%, A\, Al, and A% in the orthorhombic 
point group C2« for the cerium ion in CeRu2Alio and 
the CEF splits the J=5/2 and 7/2 multiplets of Ce 3+ 
into seven Kramer's doublets. Each of these states can 
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FIG. 1: (Color online) Measured and simulated linear polar- 
ized XAS spectra of CeRu2Alio at the M^s absorption edge at 
T = 15 K. The two simulations reproducing the experimental 
data are shown with the corresponding spatial distributions 
of the 4/ electrons as insets. The arrows indicate the spectral 
weight coming from the 4/° contribution in the initial state. 



generally be represented in the basis of \J,J Z ) as 

]T a. h \5/2,J z ) + /3jJ7/2,J z ) 

J — _2 _2 1 J — _I _J> 1 

•Jz— 2 ' 2'---'2 <Jz— 2' 2'---'2 

for values of J z which fulfill J^j ( a j + 0j ) = 1 an d 
AJ Z = ±2. The full multiplet routine takes the higher 
J =1/2 multiplet into account. However, for most Cc 
compounds the J=7/2 does not mix with the lower 
J=5/2 multiplet because the CEF splitting AEqef is 
much smaller than the spin orbit splitting AEso of 
~ 280 meV. In such a case only the A™ parameters with 
k < 4 affect the J z mixing of the lower J=5/2 multi- 
plet and it is justified to set the higher order parameters 
to zero (Af=0). For AE so > AE C ef and A^=0 the 
remaining A™ can be converted to Stevens parameters 

The bottom curves in Fig. [1] are the measured low- 
temperature linear polarized XAS spectra of CeRu2Alio. 
The main absorption features are due to the absorption 
process 3d 10 4f 1 3d 9 4/ 2 (M 4 and M 5 edges). The spec- 
tra show a clear polarization dependence in all three crys- 
tallographic directions as expected for an orthorhombic 
compound. From INS experiments the splitting between 
the CEF ground state and the first excited CEF state is 
known to be ~30meVji£ so that data taken at T = 15 K 
are representative for the CEF ground state. There is 
some additional faint spectral weight at about 888 and 
906 eV in the absorption spectra (indicated by the arrows 
in Fig. [J) due to the 3rf 10 4/° ZifiAf 1 absorption pro- 
cess which is representative for the 4/° contribution in 
the initial state. However, we note that for CeRu2Alio 
this contribution is fairly small (compared with the / 
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|7/2,±3/2} 


=F0.02(1) 


±0.02(1) 


|7/2, T 7/2) 


±0.03(1) 


±0.00(1) 



TABLE I: The J z coefficients a Jz and f3j z of the two GS wave 
functions which describe the XAS data. The wave functions 
have been calculated with the full multiplet routine, using 
c as quantization axis. The corresponding A™ are given in 
Table |IT] 
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Al 


Al 


Al 


At 


\GS#l) 


4(5) 


36(2) 


110(5) 


-30(10) 


-69(5) 


|GS#2> 


-5(2) 


38(5) 


-40(10) 


-112(5) 


0(5) 



TABLE II: Crystal-field parameters A™ in meV from full mul- 
tiplet calculations with Aq = for m — 0, 2,4, 6. For con- 
version to Stevens parameters B™ see Ref. 1381 . Note that 
Stevens parameters will not yield the J=7/2 contributions to 
the wave functions listed in Table U 



in 1:2:2 compounds),— which is in agreement with re- 
cent findings from 3c? photoemission spectroscopy^ The 
small 4/° contribution implies that the impact of the 
Kondo effect on the polarization of the 4/ 1 absorption 
edges is minor. 

The full multiplet routine yields two different wave 
functions (in the following referred to as \GS#1) and 
\GS#2)) that describe the experiment very well (see 
Fig. [I}. We use the c-axis as quantization axis and 
the corresponding J z coefficients are listed in Table Q] 
Both solutions are mainly composed of the J=5/2 mul- 
tiplet and have some small admixtures from the higher 
J=7/2 multiplet. The effect of the J=7/2-contributions 
on the LD is negligible and cannot be seen in the spectra 
when changing into Stevens approximation. For solu- 
tion \GS#1) the contribution |5/2, =p3/2) is very strong 
in contrast to solution \GS#2), which is dominated by 
|5/2,±l/2) and |5/2,±5/2). The orbitals in the in- 
sets of Fig. [1] represent the spatial distributions of the 
4/ electrons for the respective GS wave functions. It 
turns out that we can find CEF parameters for both 
wave functions which satisfy the results of INS, i.e. they 
yield the CEF transition energies AEi = 30meV and 
AE2 = 46 meV and the ratio of the inelastic neutron in- 
tensities I1/I2 — 1.35.~ The respective CEF parameters 
are summarized in Table [XXJ 

In Fig. [3] the experimental and the simulated temper- 
ature dependence of the XAS spectra are shown. Both 
CEF scenarios are based on the same energy splittings 
and at 300 K the population of the excited states at 
30 and 46 meV amounts to only 20 and 10%. Con- 
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FIG. 2: (Color online) Left panel: Measured temperature- 
dependent XAS spectra of CeRu2Alin at the M5 edge. Right 
panel: Simulation of the M5 edge temperature dependence 
for \GS=ft=l). The spectra for |GS#2) show the same behavior 
and are omitted. 



sequently both CEF models yield the same qualitative 
temperature dependence. We therefore omitted the sim- 
ulations for \GS#2) in Fig. [2j As the isotropic spectrum 
isotropic = ^_E||a±^_E||&±^_E||c is temperature-independent, 
it is sufficient to consider the temperature dependence of 
two polarizations and for reasons of clarity only the M5 
absorption edges for E\\b and E\\c are included in the 
illustration. After the first cleave we measured at 15, 
40, 100, 150, 200 and 300 K, then recleaved at 300 K and 
measured in the reversed temperature order. Moreover, 
we cross-checked and reproduced the mutual polarization 
for the two available samples. There is also no change of 
polarization up to 100 K which shows that no state gets 
populated, i.e. there is indeed no low lying crystal-field 
excitation, which is in accordance with INS resultsi&iii 
The change of polarization at 300 K with respect to 15 
K is in agreement with excited CEF states at 30 meV 
and above (see right panel in Fig. [2]). 

The combined analysis of XAS and INS data still does 
not yield a unique solution for the ground-state wave 
function. We therefore calculate the temperature depen- 
dence of the static susceptibility at H = 1 T in all three 
crystallographic directions and compare it to the experi- 
mental results found by Yutani et al^ For this purpose 
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FIG. 3: (Color online) Comparison of the experimental mag- 
netic susceptibility of CeRu2Alio as reported by Ref.0 (sym- 
bols) and the simulated ones (T > To) from the CEF param- 
eter sets belonging to \GS#1) (lines in the main panel) and 
|G5#2) (inset). 

the full Hamiltonian, including the crystal field, the mag- 
netic field and the spin-orbit coupling, is diagonalized 
and the magnetization is calculated as the Boltzmann 
weighted expectation value of the corresponding operator 
to obtain Xcef i v ~ a i b> c )- The simulation for \GS#1) 
is plotted in the main panel of Fig.j3] (lines) together 
with the data points taken from Ref. Ilo . The simulation 
for |G5#2) is included as inset. For |GS#1) both the 
anisotropy and the qualitative temperature evolution are 
nicely reproduced for T > T . The 1 5/2, J z ) contributions 
of solution |GS#1) agree quite well with the CEF analy- 
sis in Stevens approximation by Yutani et alJ£ and the- 
oretical studies by Hanzawa who determined the 4/ level 
structure of CeRu2Al 10 in a point-charge model^ Note 
that Hanzawa uses the 6-axis as quantization axis. The 
change of coordinate system from abc \\ yzx to abc \\ xyz 
can be carried out according to Ref. HH. In contrast, solu- 
tion |G5#2) does not match the measured susceptibility 
well. It yields a cross-over of Xcef an d Xcef a ^ about 
100 K so that we exclude it as a possible ground-state 
wave function. 



The CEF ground-state wave function \GS#1) yields 

the magnetic moments fJ-cp° — (1.44,0.19,0.38) ^b- Here 
an infinitesimal small temperature and magnetic field 
along the respective directions a, b, or c have been taken 
into account. These moments agree very well with val- 
ues from magnetization at 55 T and 1.3 K, which find 1.3 
and 0.3 \ib along the easy axis a and the ordering axis c, 
respectively,-- Also, the small moment of \x c CF = 0.38 [is 
agrees well with the experimental findings for the ordered 
magnetic moment,— indicating that CEF effects are 
responsible for the small size of [i or d- 

The remaining discrepancy between the measured and 
CEF-only susceptibility |G5#1) in Fig. [3] can be ex- 
plained by considering corrections due to molecular 
and/or exchange fields and anisotropic hybridization be- 
tween 4/ and conduction electrons. Such an anisotropic 
hybridization has been observed in optical conductivity 
measurements^ and the resulting anisotropic Kondo in- 
teractions will affect the anisotropy of the static suscep- 
tibility. 

In conclusion, the full multiplet simulation yields two 
GS wave functions which reproduce the low-temperature 
XAS data very well. By combining our XAS with suscep- 
tibility results^ we unambiguously identify |G5#1) (see 
Table |T]) as the GS wave function for CeRi^Alio, which 
is in good accordance with the theoretical findings by 
Hanzawa We can further give a set of CEF parameters 
that satisfies XAS, susceptibility, and INS data, 10 and the 
resulting moments are in accordance with moments from 
magnetization measurements^ The fact that the CEF- 
only moment along c amounts to only pf CF = 0.38 \ib im- 
plies that the small value of the ordered moment^— ^ 
can be explained solely with CEF effects. The absence 
of any temperature dependence in the XAS data up to 
100 K confirms that the excitation at 8 meV in the inelas- 
tic neutron spectral is not a CEF transition. The above 
results show that the selection rules for linear polarized 
light make soft XAS a powerful tool to determine the 
CEF ground state wave function of orthorhombic Kondo 
insulators even in the presence of spin gaps. 

This work was supported by DFG grant AOBJ 583872, 
Germany and KAKENHI No. 20102004 of MEXT, 
Japan. 
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